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ABSTRACT 

We present results based on Y JI\ S photometry of star clusters located in the outer¬ 
most, eastern region of the Small Magellanic Cloud (SMC). We analysed a total of 51 
catalogued clusters whose colour-magnitude diagrams (CMDs), having been cleaned 
from field-star contamination, were used to assess the clusters’ reality and estimate 
ages of the genuine systems. Based on CMD analysis, 15 catalogued clusters were found 
to be possible non-genuine aggregates. We investigated the properties of 80% of the 
catalogued clusters in this part of the SMC by enlarging our sample with previously 
obtained cluster ages, adopting a homogeneous scale for all. Their spatial distribution 
suggests that the oldest clusters, log(i yr -1 ) ^ 9.6, are in general located at greater 
distances to the galaxy’s centre than their younger counterparts - 9.0 ^ log(t yr -1 ) ^ 
9.4 - while two excesses of clusters are seen at log(t yr -1 ) ~ 9.2 and log(t yr -1 ) ~ 9.7. 
We found a trail of younger clusters which follow the Wing/Bridge components. This 
long spatial sequence does not only harbour very young clusters, log(f yr -1 ) ~ 7.3, 
but it also hosts some of intermediate ages, log(f yr -1 ) ~ 9.1. The derived cluster and 
field-star formation frequencies as a function of age are different. The most surprising 
feature is an observed excess of clusters with ages of log(t yr -1 ) < 9.0, which could 
have been induced by interactions with the LMC. 
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1 INTRODUCTION 


The near-infrared VISTA Survey of the Magellanic Clouds 
(VMC) system (Cioni ~et al.|201 1) has been designed to ob¬ 
tain three epochs of data in the near-infrared passbands Y 
and J, and 12 epochs in the K s , in order to reach a nominal 
survey depth of Y = 21.9, J = 21.4 and K s = 20.3 mag (at a 
signal-to-noise ratio of 10) for individual tiles of ~ 1.5 deg 2 


* Based on observations obtained with VISTA at ESO under pro¬ 
gramme ID 179.B-2003. 
f E-mail: andres<Q!oac.uncor.edu 


in size. The VMC survey covers a total area of ~ 170 deg 2 , 
comprising the Large Magellanic Cloud (LMC), the Small 
Magellanic Cloud (SMC), the Magellanic Bridge and a few 
tiles covering the Magellanic Stream. Individual epochs have 
exposure times of 800 s (Y and J) and 750 s ( K s ). The multi¬ 
epoch observations allow us to minimize variability effects in 
colours, particularly for bright objects, and guarantee homo¬ 
geneous observing conditions among different passbands for 
a given epoch. Seeing constraints, imposed for the purpose 
of homogenizing crowded and uncrowded field observations, 
range between 1.0” and 1.3” ( Y ), 0.9” and 1.2” (J), and 
0.8” and 1.0” ( K a ) and may exceed those values by 10% ac- 
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cording to observing policy. The averaged K s magnitudes re¬ 
sulting from the multi-epoch observations for RR Lyrae and 
Cepheid variable stars, with an accuracy of a few hundredths 
of a magnitude, will allow us to unveil the 3D structure of 


ure 1 shows the distribution of the VMC tiles superimposed 


the Magellanic system (Ripepi et al. 2012 Moretti et al. 


2014). On the other hand, stacking of these observations al¬ 


lows us to detect stars at the oldest main-sequence turn-off 
(MSTO) and derive the star-formation history (SFH) across 
the system with unprecedented quality (Rubele et al .|2012 1. 

VMC data also allow us to perform a thorough study of 
the Magellanic Cloud (MC) star cluster population based 
on a homogeneous determination of star cluster parame¬ 
ters (Piatti et al. 2014]). In addition, we can firmly establish 
whether the field star population has experienced the same 
or a similar SFH as the star cluster system. The addition of 
near-infrared photometry to existing optical photometry is 
aimed at disentangling the scatter that currently exists in 
the age-metallicity distribution (e.g. Piatti|2011b Piatti & 
|Geisler|[2013[ ), also enabling us to constrain extinction vari¬ 
ations more precisely. Wide-Held VMC data will produce a 
complete census of the cluster population that will be used 
to draw statistically robust conclusions about their differ¬ 
ences between the Clouds. In addition, the search for new 
clusters based on stellar surface density distributions will 
work better in the Clouds than in the Galaxy because of the 
lower confusion along the line of sight. 

The outermost, eastern region of the SMC and the asso¬ 
ciated onset of the Bridge are of significant scientific interest 
in the context of a putative scenario involving interactions 
between the MCs, which may be traced spatially and tem¬ 
porally. Indeed, the purpose of this paper is to present a 
photometric analysis of the catalogued clusters located in 
that region based on the most complete VMC data set to 
date. The available photometric data allow us to confirm 
the physical reality of the catalogued star clusters and (for 
the genuine clusters) estimate their fundamental parame¬ 
ters. We also investigate the cluster frequency in the region 
of interest to assess whether star clusters and field stars have 
evolved as a coupled system. 

This paper is organised as follows. The VMC data col¬ 
lection and reduction is presented in Section 2, while the 
cluster sample is described in Section 3. Sections 4 and 5 
deal with the cleaning procedure pertaining to the clusters’ 
colour-magnitude diagrams and the astrophysical property 
estimates for the clusters, respectively. In Section 6 we dis¬ 
cuss the photometric results, derive the cluster frequency 
and compare it to the field star formation history. We sum¬ 
marize our conclusions in Section 7. 


VMC DATA COLLECTION AND 
REDUCTION 


The VMC survey (Cioni et al. 2011) has collected more than 
half of the observations scheduled as part of a global observ¬ 
ing campaign with the VISTA telescope, where observations 
for five other ESO public surveys are also carried out in ser¬ 
vice mode ( |Arnaboldi et al.||2013[ ). The survey’s SMC and 
Magellanic Bridge areas will eventually cover 42 deg 2 and 
21 deg 2 (27 and 13 VISTA tiles), respectively. One tile cov¬ 
ers uniformly an area of ~ 1.5 deg 2 , representing a mosaic 
of six paw-print images in a given passband ( YJK a ). Fig¬ 


on that of the star clusters (dots) catalogued by Bica et al. 
(20081 hereinafter B08) in the SMC and the western Bridge 


regions. On the other hand, viewing the SMC as a triaxial 
galaxy (and adopting the declination, Dec, right ascension, 
RA, and line-of-sight as the three axes), Crowl et al. (2001) 


found axial ratios of approximately 1:2:4. Based on this re¬ 
sult, and to describe the clusters’ spatial distribution, Piatti 


et al. (2007a) used an elliptical rather than a spherical frame¬ 


work so as to reflect more meaningfully the flattening of the 
galaxy. This reference system - with one of its axes parallel 
to the SMC Bar and the other perpendicular to that direc¬ 
tion - seems more appropriate to describe the clusters’ age 
and metallicity distributions than one with axes along the 


right ascension and declination directions (see, e.g., Piatti 
et al.|2008 Piatti|2011b 2012b I. Two ellipses with semima¬ 


jor axes of 2.4° and 6.0° are also shown. The ellipses are 
centred at RA = O Oh 52m 45s, Dec = -72° 49' 43" (J2000) 
(Crowl et al.[20 01) and have axis ratios, b/a = 1/2. Thus the 
clusters’ spatial distribution correlates better with a pseudo¬ 
elliptical (projected) distance coordinate measured from the 
galactic centre than with the radial distance, or with dis¬ 
tances defined along the right ascension or declination axes. 

The VMC data are processed with the VISTA Data 
Flow System pipeline, version 1.3 (VDFS Irwin et al.|2004 1, 
and calibrated to match the VISTA photometric system, 
which is close to the Vegamag system. We extracted the 
observational data from the VISTA Science Archive (VSA, 
Cross et al. 2012| . The processed paw-print images were 


used to derive the effective point-s pread fun ctions ( P SFs) 
using the irai0/daophot routines (Stetson, Davis & Crab- 


tree|1990). We generated a reference PSF, which was con¬ 


volved with the paw-print images to homogenize the result¬ 
ing PSFs. We repeated these steps for each epoch separately. 
Finally, all homogenised paw-print images were combined 
using the swarp tool (Bertin et al. 20021, as described in 
|Rubele et al.| ( f2012| ) , thus generating deep tile images with 
homogeneous PSFs. 

We performed PSF photometry on the homogenized, 
deep tile images of VMC tiles SMC 3_5, 4 5, 5_6, 6_5 and BRI 
2_3. We focussed on these tiles because they are among the 
first fully or mostly completed tiles in the eastern region of 
the SMC (see column 4 of Table 1). Figure 1 shows that they 
extend over a vast region on the outermost, eastern side of 
the SMC and the associated onset of the Bridge, where early 
galactic stellar populations and those recently formed are 
supossed to co-exist as witnesses of the galaxy’s formation 
and its tidal interaction with the Large Magellanic Cloud 
(LMC). Note that these tiles include, statistically speak¬ 
ing, an important percentage of the star clusters catalogued 
in that region, thus allowing us to draw robust conclusions 
about their star-formation history. The PSF model was cre¬ 
ated using ~2500 stars which were uniformly distributed 
and had magnitudes close to the saturation limit + 1.5 mag 
(for the VMC survey, the single paw-print saturation lim¬ 
its are Y = 12.9 mag, J = 12.7 mag and K s = 11.4 mag). 


1 IRAF is distributed by the National Optical Astronomy Ob¬ 
servatories, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under contract with the U.S. 
National Science Foundation. 
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Subsequently, we used the allstar routine to perform the 
final PSF photometry on the deep tile images in all three fil¬ 
ters and correlated the resulting tables, adopting a tolerance 
of one arcsec. We applied apertu re corrections using cata- 
logues retrieved from the VSA|^ (Lewis, Irwin & Bunclark 


2010 Cross et al.|2012l for the bulk of the observed stars. 


We performed a large number of artificial-star tests to 
estimate the incompleteness and error distribution of our 
data for each deep tile image. For each VMC tile we gener¬ 
ated ~ 20 x 10 6 artificial stars following the steps outlined 


Rubele et al. (2012). We used a spatial grid of 25 pixels 


in width and with a spatial magnitude distribution propor¬ 
tional to the square of the magnitude. This latter choice 
allowed us to better map the completeness and error levels 
in the less complete regions of the colour-magnitude dia¬ 
gram (CMD). Table 1 lists the magnitudes of stars with 
photometric errors of less than 0.10 mag (mago.i) and the 
magnitudes representing the 50% completeness level. Our 
previous experience ( |Cioni et al.||20lT) |Rubele et al. |[2012[ 
Li et al.]|2014l taught us that the widest colour range, i.e., 


the Y — K s colour, is optimal for cluster studies, because 
it makes it easier to distinguish different cluster main se¬ 
quences (MSs), particularly their turn-off (TO) regions, and 
the red-giant phases. This colour is also characterized by a 
higher sensitivity to reddening and metallicity than either 
the Y — J or J — K s colours. Consequently, the subsequent 
analysis is mainly based on the K s versus Y — K s CMDs; 
the J versus Y — J and K s versus J — K s CMDshave been 
constructed and used to confirm the isochrones matched in 
the Kg versus Y — K B CMDs. 


3 THE CLUSTER SAMPLE 

We delimited the outermost eastern region of the SMC and 
the associated onset of the Bridge by an ellipse with a semi¬ 
major axis of 2.4° and relative RA coordinates ^ 1.7° (see 
Fig. 1). The star cluster population in that region com¬ 
prises 79 objects according to the B08 catalogue; 27 (34%) 
have been studied in previous photometric studies (see, e.g., 
|Glatt, Grebel fe Koch|2010||Piatti et al.|2011| ). Fifty-one out 
of the 79 catalogued clusters are spread across the VMC tiles 
considered here, which added to 12 other previously studied 
clusters located outside these VMC tile areas, results in a 
more significant cluster sample (80%) than that previously 
studied. Note that the 51 clusters on the VMC tiles do not 
only represent a larger percentage of the cluster population 
in that region (65%) with respect to the previously studied 
clusters, but here we also analyse them in a homogeneous 
manner, thus enabling us to trace their formation history 
spatially. To date, the catalogue of B08 is the most com¬ 
plete compilation of star clusters in the SMC covering the 
extension of both Clouds and the Bridge. Particularly, be¬ 
cause of the lower surface brightness of the background and 
the less-populous nature of the outer SMC regions and the 
Bridge, star clusters stand out more, and incompleteness ef¬ 
fects are expected to be less important than in the inner 
regions of the Clouds. 

Recognising catalogued star clusters in deep VMC tile 


2 http://horus.roe.ac.uk/vsa/ 


images is neither straighforward nor simple. On the one 
hand, the catalogued objects were originally identified from 
optical images (e.g., from the Digitized Sky Survey, DS^J 
images) which sometimes look rather different compared 
with their appearance at near-infrared wavelengths. In addi¬ 
tion, the spatial resolution and depth of the images on which 
the clusters were identified differ from the equivalent pa¬ 
rameters pertaining to the VMC images. Thus, for instance, 
single relatively bright stars might look like an unresolved 
compact cluster in images of lower spatial resolution, or un¬ 
resolved background galaxies could be mistaken for small 
star clusters in shallower images. Offsets in the coordinates 
compiled by B08 with respect to the objects’ centres cannot 
be ruled out either. 

To avoid mismatches between observed objects and the 
actual list of catalogued clusters, we first overplotted the 
positions of the catalogued clusters (B08) on the deepest K s 
image. This way, based on using the coordinates provided 
by B08, we visually recognized the observed clusters one by 
one in the K s image. Next, we searched for these clusters 
in the DSS and downloaded 15' xl5' B images centred on 
the coordinates resolved by the SIMBAL^astronomical data 
base and compared them with equivalent cut-outs derived 
from the K s VMC survey. Thus, we correctly recognized the 
clusters in both the optical and near-infrared regimes. Note 
that the main aim of this task is to confirm the cluster co¬ 
ordinates and sizes, in order to extract from the tile-image 
PSF photometry the magnitudes of the stars in the clus¬ 
ter region. We are not interested in properties such as the 
clusters’ structure, stellar density profiles or radii, but in 
the stars which allow us to meaningfully define the clusters’ 
fiducial sequences in their CMDs. Table 2 lists the complete 
cluster sample, as well as the coordinates and radii adopted 
for extracting the stellar PSF photometry; the cluster’s ra¬ 
dius was taken either from a visual inspection of the deepest 
K s image (where the profile disappears into the background 
noise), from B08 or from both sources combined. The ob¬ 
served objects are of small angular size, typically ~ 0.7 / in 
diameter (~ 12.2 pc). 


4 CLEANING THE CLUSTER CMDS 

In general, the extracted cluster CMDs represent the com¬ 
posite stellar populations distributed along the respective 
lines of sight. Since they account for the luminosity func¬ 
tion, colour distribution and stellar density towards a given 
region on the sky, CMD analysis alone might lead to an in¬ 
correct interpretation. On the other hand, the CMDs of the 
stars located within a region around the catalogued cluster 
centres are a helpful tool to assess the reality of the density 
enhancements. They may imply that we are dealing with 
the presence of a genuine star cluster, a chance grouping of 


3 The Digitized Sky Surveys were produced at the Space Tele¬ 
scope Science Institute under U.S. Government grant NAG W- 
2166. The images of these surveys are based on photographic data 
obtained using the Oschin Schmidt Telescope on Palomar Moun¬ 
tain and the UK Schmidt Telescope. The plates were processed 
into the present, compressed digital form with the permission of 
these institutions. 

4 http://simbad.u-strasbg.fr/simbad/ 
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stars along the line of sight or a non-uniform distribution 
of interstellar material in the region of interest. Note that 
catalogued cluster candidates appear on the sky as small 
concentrations of stars, although this does not necessarily 
imply that such concentrations constitute real, physically 
bound systems. 

For these reasons, we first statistically constructed 
CMDs representing the field along the line of sight towards 
the individual clusters, which we then used to clean the clus¬ 
ter CMDs. We employed the cleaning procedure developed 


by Piatti & Bica (2012 see their fig. 12). The method com¬ 


pares the extracted cluster CMD to four distinct CMDs com¬ 
posed of stars located reasonably far from the object, but 
not too far so as to risk losing the local field-star signa¬ 
ture in terms of stellar density, luminosity function and/or 
colour distribution. The four field regions were designed to 
cover an equal area as that of the cluster (a circular area of 
three times the cluster radius), and they were placed to the 
north, east, south and west of the cluster area. Note that 
large areas were chosen to increase the statistics and hence 
the reliability of the cleaning process. 

Comparisons of field and cluster CMDs have long been 
done by comparing the numbers of stars counted in boxes 
distributed in a similar manner throughout both CMDs. 
However, since some parts of the CMD are more densely 
populated than others, counting the numbers of stars within 
boxes of a fixed size is not universally efficient. For instance, 
to deal with stochastic effects at relatively bright magni¬ 
tudes (e.g., fluctuations in the numbers of bright stars), 
larger boxes are required, while populous CMD regions can 
be characterized using smaller boxes. Thus, use of boxes of 
different sizes distributed in the same manner throughout 
both CMDs leads to a more meaningful comparison of the 
numbers of stars in different CMD regions. By starting with 
reasonably large boxes - (AA' S , A(Y — K s )) = (1.00, 0.25) 
mag - centred on each star in the four held CMDs and by 
subsequently reducing their sizes until they reach the stars 
closest in magnitude and colour, separately, we defined boxes 
which result in use of larger areas in held CMD regions con¬ 
taining a small number of stars, and vice versa (see Fig. 2, 
top right-hand panel, where we used an annulus -outer and 
inner radii equal to 3.163 and 3.0 times the cluster radius- 
around the cluster instead of one of the four selected cir¬ 
cular areas for illustrative purposes). Next, we plotted all 
these boxes for each held CMD on the cluster CMD and 
subtracted the star located closest to each box centre. 

Since we repeated this task for each of the four held- 
CMD box samples, we could assign a membership proba¬ 
bility to each star in the cluster CMD. This was done by 
counting the number of times a star remained unsubtracted 
in the four cleaned cluster CMDs and by subsequently divid¬ 
ing this number by four. Thus, we distinguished held popu¬ 
lations projected onto the cluster area, i.e., those stars with 
a probability P ^ 25%, stars that could equally likely be as¬ 
sociated with either the held or the object of interest (P = 
50%), and stars that are predominantly found in the cleaned 
cluster CMDs (P ^ 75%) rather than in the held-star CMDs. 
Figure 2 illustrates the performance of the cleaning proce¬ 
dure for the HW77 held, where we plotted three different 
CMDs: that for the stars located within the cluster radius 
(top left-hand panel), a single-held CMD for an annulus - 
outer and inner radii equal to 3.163 and 3.0 times the cluster 


radius- centred on the cluster, as well as the cleaned cluster 
CMD (bottom left-hand panel). The schematic diagram with 
a superimposed circle of radius equal to the cluster radius 
is shown in the bottom right-hand panel. The pink, light 
and dark blue solid circles in the bottom panels represent 
stars with P ^ 25%, P = 50% and P ^ 75%, respectively. 
Note that when comparing observed cluster and held CMDs 
for the HW 77 region (top left- and right-hand panels), the 
cluster region contains brighter red clump (RC) stars and a 
younger MSTO : A' s (RChw 77 ) ~ 17.0 mag, AT s (RC fle id) ~ 
17.5 mag, A' s (MSTOhw77) ~ 19.5 mag and A' s (MSTOfi e id) 
~ 20.5 mag, which suggests that HW77 is projected against 
a relatively older composite stellar population. A full set of 
hgures for the remaining objects studied here is provided in 
the online version of the journal. 

To assess whether a catalogued cluster is a genuine ag¬ 
gregate, we inspected both the cluster’s schematic diagram 
and the cleaned CMD. Objects showing some apparent con¬ 
centration in the sky of stars with P ^ 75% that do not 
dehne clear cluster sequences in the cleaned CMDs were as¬ 
sumed to be possible non-genuine aggregates. We found that 
nearly 30% (15 catalogued clusters) of our cluster sample 
fall into this latter category. This relatively high percentage 
compares well with recent results of different cluster sample 
analyses in both Magellanic Clouds, which suggest that the 
B08 cluster catalogue might contain between 10% and 30% 
of possible non-genuine aggregates (Piatti fc Bica|2012 Pi 


|atti||2014a| |Piatti et al.| [2014). Likewise, we found two clus¬ 
ters (GKH22 and GKH24) that could not be resolved by 
our photometry and another (BS226) that could not be rec¬ 
ognized in any image using the coordinates given by B08. 
The relevant information about the status of these objects 
is provided in the final column of Table 2. 


5 FUNDAMENTAL CLUSTER PARAMETERS 

We used the cleaned cluster CMDs to estimate the fun¬ 
damental cluster parameters by matching the observations 
with the theoretical isochrones of Bressan et al.| ( |2012[ ). In 
performing this task, one has to deal with parameters such 
as the reddening, distance, age and metallicity. Our strategy 
consisted of obtaining the reddening values from an indepen¬ 
dent source, assuming the mean SMC distance modulus for 
all clusters and adopting for the cluster ages the ages of 
the isochrones which best reproduced their CMD features. 
We started with isochrones for a metallicity of Z = 0.003 
([Fe/H] = —0.7 dex), which corresponds to the mean SMC 
cluster metal content during the last ~ 2-3 Gyr ( |Piatti fc| 
Gcislcr 2013), and employed isochrones for other metallici- 
ties as required. Note that the difference in Y — K s colour 
between isochrones for Z = 0.003 and 0.001 ([Fe/H]= —1.1 
dex) is ~ 0.1 mag in the MSTO region, ~ 0.05 mag at the 
subgiant branch (SGB) and ~ 0.15 mag at the tip of the 
red-giant branch for ages between log(i yr -1 ) = 9.3 (2 Gyr) 
and log(t yr -1 ) = 9.6 (4.0 Gyr). This difference is ~ 0.05 
mag smaller for log(t yr -1 ) = 9.8 (6.3 Gyr). This small dif¬ 
ference suggests that the Y — K s colour is not sensitive to 
metallicity differences smaller than A[Fe/H]=0.4 dex if we 
keep in mind the relatively sparse nature of the majority 
of our clusters and the intrinsic spread of the stars in the 
CMDs. Indeed, we tried using isochrones with metallicities 
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[Fe/H] = -0.5 and -0.9 dex and found negligible differences 
with respect to that of [Fe/H] = —0.7 dex. Nevertheless, only 
four clusters in our sample (L109, L110, B168 and HW85) 
appeared older than log(t yr - ) = 9.3, and we carefully paid 
attention to the effects of metallicity differences. 

According to Dutra et al. (2001), the SMC is opti¬ 
cally thin, characterized by average foreground and internal 
E(B — V) colour excesses of 0.05±0.05 mag and 0.04 mag 
(computed from a reddening for the SMC surroundings of 
E{B — V) = 0.02±0.02 mag), respectively. We took advan¬ 
tage of the Magellanic Cloud extinction values based on RC 
stellar photometry provided by the OGLE III collaboration 


(Udalski 2003), as described in Hasclike, Grebel & Duffau 


(2011), to estimate E(B — V)ogle in colour excesses for 18 
clusters in our sample which are located in the area covered 
by OGLE III. Likewise, we made use of the NASA/IPAC 
Extragalactic Data bas^j (NED) to infer Galactic fore¬ 
ground reddening values, E(B — V)ned, for the same clus¬ 
ter list. Next, we computed the difference, A E(B — V) = 
E(B — V)ogle hi — E(B — V)ned = 0.025 ± 0.005 mag, in 
good agreement with the mean internal reddening given by 


Dutra et al. (2001). Note that this mean SMC internal red¬ 


dening renders the isochrones 0.02 mag redder in Y — K a and 
0.01 mag fainter in K s , respectively, in the cleaned cluster 
CMDs. In matching isochrones, we started by adopting ei¬ 
ther E(B — V)ogle hi values or E(B —V)ned + 0.03 mag, 
combined with the equations E(Y — K s ) = 0.84 x E(B — V) 
and K s = Mk b + (m — M)q + 0.372 x E(B — V), assum¬ 
ing Ry = 3.1 (Cardelli, Clayton fe Mathis||1989 |Gao et ah 


2013). The final E(B — V) values are listed in Table 2. We 


found two clusters which may be affected by differential 


reddening, i.e., HW74 (Bratsolis, Kontizas & Bellas-Velidis 


2004) and HW81 (embedded in an Hu region Testor et al. 


2010). Two other clusters (BS189, BS190) are characterized 


by noticeably large E(B — V) values, that may also be in¬ 
fluenced by differential reddening. In particular, BS189 is 
by far the most reddened cluster known in the Magellanic 
Clouds; it is even more heavily reddened than any cluster in 


the highly reddened LMC/30 Doradus region (Tatton et al. 
|M3l ). 

For all clusters, we adopted the mean SMC di stance 
modulus (m—M)o = 18.90 ± 0.10 mag (GO.O/^'g kpc) (Glatt, 


Grebel & Koch 


2010) and an average line-of-sight SMC disc 


depth of 6 kpc (Crowl et al.|200T|). Bearing in mind that any 


cluster in the present sample could be placed in front of or 
behind the SMC, we conclude that the difference in apparent 
distance moduli could be as large as A (K s —Mk s ) ~ 0.2 mag. 
This difference is much smaller than the difference between 
two closely spaced isochrones as used here (A log(t yr -1 
= 0.1, A Mk b ~ 0.4 mag), so that adoption of a unique 
value for the distance modulus does not dominate the final 
error budget incurred in matching isochrones to the cluster 
CMDs. 

In the matching procedure, we used subsets of 
isochrones ranging from A log(t yr -1 ) = —0.3 to +0.3, strad¬ 
dling the initial rough age estimates. Since the |Bressan et ak| 
12012) isochrones are defined in the Vegamag system (where 


5 http://ned.ipac.caltech.edu/. NED is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under 
contract with NASA. 


Vega has a magnitude of zero in all filters, by definition), 
we subtracted 0.074 mag in Y and 0.003 mag in I\ B from 
the isochrones before matching them to the cluster CMDs 
(Rubele et al. 20151. Thus, we worked with both theoreti¬ 
cal models and observed CMDs which were defined in the 
VISTA system. We found that isochrones bracketing the de¬ 
rived mean age by A log(t yr -1 ) = ±0.1 represent the over¬ 
all age uncertainties owing to the observed dispersion in the 
cluster CMDs. Although in some cases the age dispersion 
is smaller than A log(t yr - ) = 0.1, we prefer to keep the 


former value as an upper limit to our error budget (Piatti 
et al.|2011 Piatti|[2014b[ Piatti et al.||2014| among others). 


Finally, we adopted for the cluster age/metallicity those val¬ 
ues corresponding to the isochrone which best reproduced 
the cluster’s main features in the CMD. Table 2 lists the 
resulting age and metallicity values for 33 of the confirmed 
clusters, while the bottom left-hand panel of Fig. 2 illus¬ 
trates our isochrone-matching procedure in the context of 
the CMD of HW77. 

Fourteen clusters have previous age estimates based on 
matching of theoretical isochrones to their CMDs. We pre¬ 
viously studied seven of these based on Washington CTi 
photometry and the remaining clusters were analysed by 


Glatt, Grebel & Koch (20101 as part of the Magellanic Cloud 
Photometric Surveys (MCPS [Zaritsky et al.|2002[ ). We took 
advantage of the availability of the Washington data to per¬ 
form a sound comparison of the ages resulting from both the 
Washington photometry and the present results. First, we 
cleaned the cluster CMDs following the precepts described 
in Section 4. Then, we adopted the mean SMC distance 
modulus, E(B — V ) colour excesses and metallicities dis¬ 
cussed above. Finally, we matched isochrones to the cluster 
CMDs, shifted by E(C - 7\) = 1.97 E{B - V ) and Ti = 
MT 1 +2.62E(B — V) + (m—M)o ( |Geisler fc Saraj edini 1999| ). 
The resulting ages are listed in the final column of Table 2. 
When comparing these values with the VMC cluster ages, we 
found a difference of A log(i yr -1 ) = — 0.05±0.10, with the 
VMC ages being slighlty older. On the other hand, a direct 
comparison of the age estimates derived by|Glatt, Grebel fe| 


Koch (20101 (see the final column of Table 2) and our values 
shows that our ages are A log(t yr - ) = 0.5 ± 0.5 older. We 
recall that Glatt et al. did not perform any decontamina¬ 
tion of field stars from the cluster CMDs and they therefore 
had to decide whether to include luminous supergiants in 
their fits, since the apparent MSTO of sparse, young clus¬ 
ters is subject to pronounced statistical fluctuations. It is 
possible that the lack of cleaned CMDs and their shallower 
photometry did not allow them to achieve more reliable age 
estimates. 


6 ANALYSIS 

In this section we address the cluster-formation history on 
the outermost, eastern side of the SMC and the connected 
onset of the Bridge by analysing the spatial dependence of 
their age distribution. We first enlarged the cluster sample 
by adding previously studied objects located in this region 
but beyond the VMC tile areas. We found a total of 12 clus¬ 
ters with available CT\ Washington photometry and ages 
estimated from the matching of isochrones to their CMDs. 
As shown above, these age estimates are in very good agree- 
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ment with those of the VMC clusters. We needed to correct 
the CTli-based ages by only A log(t yr -1 ) = —0.05 to es¬ 
tablish a homogeneous age scale, i.e., that of the present 
VMC cluster ages. The properties of the additional cluster 
sample are listed in Table 3. 

The enlarged cluster sample is shown in Fig. 3 using 
solid circles. We used a colour scale to distinguish among 
clusters with different log(t yr -1 ) values, so that darker 
colours correspond to older clusters. The spatial distribu¬ 
tion of the clusters seems to suggest that while the 8 oldest 
clusters - log(f yr -1 ) ^ 9.6 (see Tables 2 and 3) - are in 
general located at greater distances to the galaxy’s centre 
than the younger clusters - 9.0 ^ log(f yr -1 ) ^ 9.4 - there 
is also a trail of young clusters, log(f yr -1 ) < 8.4, along the 
onset of the Bridge. Since we are dealing with 80% of the 
catalogued clusters in this part of the SMC - the remaining 
20% are mostly located in the Bridge - our results can be 
considered as based on an unbiased cluster sample. 

Viewing the SMC as a triaxial galaxy with the Dec, RA 
and line of sight as the three axes, |Crowl et al. ([2001) found 
axial ratios of approximately 1:2:4. Based on this result, 
and aiming at estimating the clusters’ projected distances 
to the SMC’s centre, we adopted the elliptical framework 
described in Section 2, where a is the semimajor axis of an 
ellipse centred on the SMC centre, for a b/a axis ratio of 1/2 
and one point of its trajectory coinciding with the position 
of the cluster of interest. For our sample clusters with log(f 
yr - ) is 9.6 we obtained a = (6.0 ± 1.3)°, in excellent agree¬ 
ment with the mean value for the 16 known SMC clusters 


intermediate-age stars that could have been stripped from 
the SMC, which has also been predicted by models of gas 


spanning the same age range i 

Dolphin et al. 

2001 

Piatti 

2011a, Da Costa & Hatzidimitriou 1998||Mighe 

1, Sarajedini 

Sz French 1998| Piatti et al.|2007c Piatti|2012a). Note that 


both the present cluster sample and the handful of the old¬ 
est known clusters span a semimajor axis range between ~ 
2° and 9°; the latter are found at any position angle. For 
clusters with 9.0 ^ log(t yr -1 ) ^ 9.4, we obtained a = (4.9 
± 1.8)°. Although the latter result is statistically indistin¬ 
guishable from that for the oldest clusters in the sample, the 
smaller mean value suggests that the intermediate-age clus¬ 
ters are spatially slightly more concentrated than the older 

Carrera| 


cluster population. This trend is in agreement with 


et al. (20081, Piatti ( 2012b|a I and Cignoni et al. 


(2013), 


who showed that SMC clusters and field stars follow a sim¬ 
ilar trend in the sense that the larger the semimajor axis is, 
the older their ages are, with a non-negligible dispersion. 

The other interesting feature of Fig. 3 is the trail of 
younger clusters that seems to follow the Wing/Bridge com¬ 
ponents. Such a long spatial sequence, which seems to orig¬ 
inate directly in the inner SMC disc, does not only harbour 


replenishment from the SMC to the LMC (Bekki & Chiba 


20071. Harris (20071 did not detect intermediate-age stel¬ 


lar populations associated with the Bridge, which led him 
to conclude that the material which was stripped from the 
SMC to form the Bridge may have consisted of very nearly 
pure gas. However, the present cluster ages agree well with 
results from field-star studies which showed that the Bridge 
does not only contain gas from which new generations of 
stars might be formed but also older stellar populations. 
Whether (i) this star-formation process took place along 
the Wing/Bridge components, (ii) clusters/field stars were 
stripped from the SMC in the direction of the LMC and 
(iii) their formation history is a continuous or a burst-like 
chain of events are questions that are still being debated and 
which deserve further investigation. From the perspective of 
the clusters’ ages, the Wing/Bridge structures seem to have 
existed at least during the past 1 -2 Gyr, i.e., log(t yr -1 ) ~ 
9.0-9.3. 

Figure 4 shows a time series pertaining to the clusters’ 
spatial distribution in intervals of Alog(t yr -1 ) = 0.2; the 
lower limit included in the interval. We have drawn, in ev¬ 
ery panel, all catalogued clusters (dots) and the same two 
ellipses as in Fig. 3. The present cluster sample (Tables 2 
and 3) is shown as solid circles. Figure 4 reveals that the 
ages of the clusters in the outermost, eastern part of the 
SMC cover the entire age range, with the exception of the 
oldest ages and the interval 7.4 ^ log(t yr - ) < 7.8. The 
lack of SMC clusters that are as old as the Galactic glob¬ 
ular clusters has long been known ( |Glatt et al. 20081. On 
the other hand, the absence of clusters in the younger age 
range indicated can be li nked to the minimum age found by 


Glatt, Grebel & Koch (2010), log(t yr ) ~ 7.5, based on an 


age distribution composed of 821 clusters and associations 
distributed across the main body of the galaxy; only ages 
with Alog(f yr -1 ) < 0.5 were used. 

Clusters older than log(t yr -1 ) ~ 9.2 appear to be found 
predominantly in the outer disc, whereas younger clusters 
tend to be mostly located along the Wing/Bridge interface. 
Two peaks populated by older clusters are seen at log(f yr -1 ) 
~ 9.0-9.4 and log(t yr - ) ~ 9.6-9.8, in very good agreement 
with enhanced episodes of cluster formation that occurred 
throughout the galaxy Piatti ( 201 lb| and reference therein). 
Both peaks were found by jNoe l et al. ([2009) from the recov¬ 
ery of the formation history of stars observed in 12 fields 
in the B, R bands with the 100 inch telescope at Las Cam- 
panas Observatory. |Harris & Zaritsky| ([2004]) also detected 
the younger peak based on their MCPS data, whereas |Rezaei| 
kh. et al. (20141 showed the presence of the older peak based 


very young clusters, log(t yr x ) ~ 7.3, but also some of in¬ 
termediate ages, log(t yr -1 ) ~ 9.1. |McCumber, Garnett fe| on observations of long-period variable stars. 


Dufourj (2005) imaged a star field located in the Wing with 
the Hubble Space Telescope (HST)’s WFPC2 camera and 
also found a wide age range, i.e., log(f yr -1 ) ~ 8.0-9.2, while 
|Sabbi et al.| ( |2009| ) used the HST /ACS camera to observe 
two fields in the Wing. They found evidence of enhance¬ 
ment in the star-formation activity starting from log(t yr -1 ) 
~ 8.7. Continuous star-formation activity at intermediate 
ages, as well as a very recent star-formation phase, have re¬ 


cently been found by Ripepi et al. (20141 and Skowron et al. 
(|2014[). In addition, pagheri, Cioni fe Napiwotzki| (|2013|) 


and Noel et al. (2013) showed evidence of the presence of 


6.1 Star cluster frequency 

We constructed the star cluster formation frequency as a 
function of age (CF) - the number of clusters per unit time 
interval as a function of age - to compare it with that de¬ 
rived under the assumption that they are characterized by a 
similar formation history as that of the field stars projected 
along their respective lines of sight. We were interested in 
obtaining the CFs for clusters located in tiles SMC 3 5, 4_5, 
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5_6 and 6_5 for which Rubele et al. (2015) obtained star- 


formation rates using the same VMC photometric data set. 

The first step in generating the CF consisted of con¬ 
sidering the age errors. Indeed, by taking into account such 
errors, the interpretation of the resulting CF can differ ap¬ 
preciably from that obtained using only the measured ages 
without accounting for their uncertainties. However, treat¬ 
ment of the age errors in the CF is not straightforward. 
Moreover, even if errors did not play an important role, bin¬ 
ning of the age ranges could also bias the results. At first 
glance, a fixed age-bin size may not be appropriate, since 
the resulting age distribution depends on the adopted age 
interval, and the age errors are typically a strong function 
of age. On the contrary, adoption of age-bin widths on the 
order of the clusters’ age errors in the interval of interest ap¬ 
pears more meaningful. This would lead to the selection of 
narrow bins (in linear age) for young clusters and relatively 
broader age bins at older ages (Piatti fc Geisler|2013 Piatti 


2014a I. To account for the effects of the age uncertainties in 


the CF while taking into account that any individual point 
in the CF may fall either in the respective age bin or in any 
of the adjacent bins, we followed Piatti (2014b). 


For our purposes, we first considered the cluster age 
range, split into bins of A log(f yr -1 ) = 0.10. On the other 
hand, each age data point and its associated uncertainty 
(cr(f)) covers a segment with a size of 2xcr(t) (which may or 
may not fall fully in one of the age bins) and has dimensions 
smaller, similar to or larger than the age bin in which it is 
found. For this reason, we weighted the contribution of each 
age data point to each of the associated age bins, with the 
sum total of all weights being equal to unity. The assigned 
weight was computed as the fraction of the age segment 
that falls within the relevant age bin. In practice, for each 
age interval we looked for clusters with ages that fall inside 
the age bin considered, as well as clusters where t ± a(t) 
could cause them to fall in the same age bin. Note that if 
t ± cr(f) causes a cluster’s age range to extend beyond the 
main age bin considered (e.g., from one bin to the next), 
then we consider that that cluster may have an age that also 
places it inside the extended age interval considered. Figure 
5 shows the resulting CF (thick solid line). This CF has 
been normalized to the total number of clusters. Note the 
discontinuity at log(t yr -1 ) = 7.4-7.8, which may be owing 
to the small number of clusters in our sample. The younger 
part of the CF is mainly based on clusters located in the 
inner tiles, SMC 3_5 and 4_5, while the remaining outer tiles 
(5_6 and SMC 6_5) mostly contribute to the older regime. 

To place the CF resulting from our small cluster sample 
in the overall context of the SMC’s cluster population, we 
also included the CF constructed on the basis of the clus¬ 
ter age determinations of de Grijs & Goodwin (2008; thick 
dashed line). The latter authors obtained best-fitting ages 
for more than 300 SMC clusters based on broad-band photo¬ 
metric measurements at optical wavelengths. They validated 
their resulting age estimates with respect to both those ob¬ 
tained by other teams using similar approaches and by refer¬ 
ence to their extensive validation of the modelling approach 
used. Despite the differences in both cluster numbers and 
photometric completeness limits between both studies, the 
overall trends seen in Fig. 5 are similar. 

The effects of small-number statistics are apparent in 
the stochastic fluctuations in the CF based on the VMC 


data set. The main difference between the CF constructed 
based on the VMC data set with respect to that resulting 
from the de Grijs & Goodwin (2008) data is the relatively 
larger number of older clusters in the VMC data. A careful 
comparison between both data sets reveals that the com¬ 
pleteness limit in terms of cluster masses of the VMC data 
sets in at much lower masses than that pertaining to the 
broad-band data set. The older, lower-mass clusters seen in 
the near-infrared VMC data had already faded to below the 
detection limit in the bluest filter used by de Grijs & Anders 
(2008). All other apparent differences between both obser¬ 
vational CFs can be traced back to small-number statistics. 
Note that the similarity between both CFs, which are nev¬ 
ertheless based on very different lower-mass limits, implies 
indirectly that the cluster mass function in the SMC is well 
represented by a power-law function down to the lowest ac¬ 
cessible masses (cf. de Grijs & Goodwin 2008). 

Finally, we compare the present, VMC-based CF to that 
obtained from the star-formation rates derived by |Rubele| 


et al. (2015) from the same VMC tiles and PSF photometry 


data set used here. They subdivided each tile into 12 sub- 
regions of 21.0'x21.5' (~ 0.12 deg 2 ) and derived their star- 
formation histories as described in detail in IRubele et al.l 


12012). Briefly, for a range of distance moduli and visual ex¬ 


tinction estimates, ‘partial models’ are derived for the entire 
range of metallicities and ages of relevance. Partial models 
are synthetic stellar populations, each covering small bins 
in age and metallicity, shifted to the desired distance and 
extinction. Each partial model is ‘degraded’ to the condi¬ 
tions of the actual observations by convolving them with 
the distributions of photometric errors and completeness. 
The linear combination of partial models that optimally 
matches the observed Hess diagrams is found by means of 


the StarFISH optimisation code (Harris & Zaritsky 20011. 


The coefficients of this linear combination of partial models 
(including the best-fitting distance modulus and extinction) 
are directly converted into the star-formation rates. 

We assumed that clusters are formed with a power- 
law mass distribution characterized by a slope of a = —2 
( Baumgardt et ah||2013 Piatti|[2014b I and with a rate that 
is proportional to the field-star formation rate determined by 


Rubele et al. (2015) for the individual subfields pertaining 


to our clusters, taking into account the corresponding uncer¬ 
tainties. We used cluster masses from log(M c i [Mg]) = 2.2 to 
log(M c i [Mg]) = 5.0 ( de Grijs fe Goodwin|2008||Glatt et al. 


2011) and normalized the resulting CF by the total number 


of clusters used, so that it can be compared directly to the 
observed distribution. Figure 5 shows the resulting model 
CF, where the uncertainties are indicated using thin solid 
lines. As can be seen, the observed and modeled CFs are 
different. For ages in excess of log(f yr -1 ) ~ 9.4 the model 
CF is higher than the observed frequency, possibly because 
we did not take into account any cluster dissolution. Ac¬ 
cording to Lada & Lada (20031, most - if not all - stars 


form in some sort of cluster. This implies that field stars are 
the result of cluster dissolution and do not originate from 
an independent formation mechanism, so that the difference 
between the observed and modeled CFs would then be the 
result of the prevailing cluster dissolution rate. 

The decrease of the star-formation activity that started 


at log(t yr 


9.0 and was interrupted by a burst-like for¬ 


mation event which occurred at log(t yr ) ~ 8.5 was also 
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documented by Harris & Zaritsky (20041, Noel et al. (20091, 


Sabbi et al. (2009) and Rubele et al. (2015). However, the 


observed CF is significantly higher for ages younger than 
log(t yr - ) ~ 9.0. Even though the model CF requires addi¬ 
tional refinements, the clear cluster excess deserves further 
analysis. If we now keep in mind |de Grijs, Goodwin fe An-| 
ders (20131’s results, who showed that there is no evidence of 


significant destruction other than that expected from stellar 
dynamics and evolution for simple stellar population models 
for ages up to 1 Gyr - log(f yr -1 ) = 9 - we speculate that 
such a cluster excess could be evidence of enhanced cluster 
formation in this part of the SMC with respect to that on the 
western side. Indeed, the number of clusters on that latter 
side is significantly smaller (see Fig. 3). Assuming that there 
probably has been no significant mixing of clusters among 
different regions, the observed excess of clusters reveals that 
the cluster formation history in the outermost, eastern re¬ 
gion of the SMC could have been influenced at some level 


by the galaxy’s interaction with the LMC. 

Such an interac- 

tion (Nidever et al.|2013 Rezaei kli. et al. 

2014) could have 

resulted in either by stripping mechanics 1 

Bagheri, Cioni & 

Napiwotzki 2013 Noel et al. 20131, or in 

situ star forma- 

tion processes (McCumber, Garnett & Dufour .2005 Ripepi 


et al. 2014), or a combination of both that caused the ob¬ 


served and modeled CFs to look different for ages younger 
than 1 Gyr. 


7 CONCLUSIONS 

In this paper, we have analysed CMDs of catalogued star 
clusters located in the outermost, eastern region of the SMC 
based on a YJI\ S photometric data set obtained by the 
VISTA VMC collaboration. We focussed on tiles SMC 3_5, 
4_5, 5 6, 6_5, and BRI 2_3, because they are among the first 
fully or mostly completed tiles in the eastern region of the 
SMC, which may also show evidence of the SMC’s interac¬ 
tion history with the LMC. We obtained PSF photometry 
for stars in and projected towards 51 catalogued clusters of 
small angular size, typically ~ 12.2 pc in diameter, which 
represents a meaningful sample size for scientific studies. 

We applied a field-star subtraction procedure to statisti¬ 
cally clean the cluster CMDs from field-star contamination 
to disentangle cluster features from those associated with 
their surrounding fields. The technique we employed makes 
use of variable cells to reproduce the field CMDs as closely 
as possible. We found that nearly 30% (15 catalogued clus¬ 
ters) of the cluster sample investigated do not have cluster 
sequences in the cleaned CMDs an were therefore assumed 
to be possible non-genuine aggregates. 

Based on matching theoretical isochrones in the VISTA 
system to the cleaned cluster CMDs, we obtained estimates 
of the cluster ages, assuming as initial guess for the clus¬ 
ter metalicity the value Z = 0.003. We took into account 
the SMC’s distance modulus as well as the individual clus¬ 
ter colour excesses. The resulting cluster ages span the age 
range 7.0 ^ log(t yr - ) ^ 9.8. This cluster sample forms 
part of the cluster data base which will result from the VMC 
survey and which will be used to self-consistently study the 
overall cluster formation history of the Magellanic system. 
Fourteen clusters have previous age estimates based on the 
matching of isochrones to their CMDs. We studied seven of 


these based on Washington CTj photometry, whereas the 
remaining clusters were studied by |Glatt, Grebel fe Koch| 
(2010). Although the present cluster ages are in very good 
agreement with those resulting from analysis of Washington 
photometry, Glatt, Grebel & Koch 20101’s ages are much 
younger, possibly because of the lack of cleaned CMDs and 
the shallower photometry used by these authors. 

We complemented the properties of 80% of the cata¬ 
logued clusters in this part of the SMC by increasing our 
sample by addition of previously determined cluster ages, 
taking care to adopt a homogeneous scale. Their spatial dis¬ 
tribution seems to suggest that the oldest clusters - log(t 
yr -1 ) ^ 9.6 - are in general located at greater distances from 
the galaxy’s centre than their younger counterparts, 9.0 ^ 
log(t yr -1 ) ^ 9.4, which is in agreement with previous stud¬ 
ies of star clusters and field stars. The latter showed that the 
greater their distances to the SMC centre, the older the clus¬ 
ters’ ages are, with a non-negligible dispersion. Within the 
older cluster population, two excesses of clusters are seen, at 
log(t yr -1 ) ~ 9.0-9.4 and log(t yr -1 ) ~ 9.6-9.8, in very good 
agreement with the timing of enhanced episodes of cluster 
formation previously found throughout the galaxy. 

We also found a trail of younger clusters that seems 
to follow the Wing/Bridge components. Such a long spatial 
sequence, which extends from the inner SMC disc, does not 
only harbour very young clusters - log(t yr -1 ) ~ 7.3 - but 
it also includes some of intermediate ages, log(t yr -1 ) ~ 9.1. 
Therefore, the Wing/Bridge structures seem to have existed 
at least during the last 1-2 Gyr, log(t yr - ). 

The observed and modeled CF, excluding tile BRI 2_3, 
are different. For ages in excess of log(t yr -1 ) ~ 9.4, the 
model CF, obtained by assuming a cluster formation rate 
proportional to that of the field stars and a power-law mass 
distribution with a slope of a = —2, is higher than the 
observed frequency, possibly because we did not take into 
account any cluster dissolution. For younger ages, the ob¬ 
served CF is - surprisingly - significantly higher than the 
model distribution, which may be partially owing to the fact 
that field-star formation activity is known to decrease from 
log(t yr -1 ) ~ 9.0 until a burst-like formation event occurred 
at log(t yU 1 ) ~ 8.5. Nevertheless, even though the model 
CF requires additional refinements, if we assume that there 
probably has been no significant mixing of clusters among 
different regions, the observed excess of clusters reveals that 
the cluster formation history in the outermost, eastern re¬ 
gion of the SMC could have been influenced at some level 
by the galaxy’s interaction history with the LMC. 
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Table 1 . VMC tile information. 


Tile ID 

RA (J2000) 
(°) 

Dec (J2000) 
(°) 

Completion in 
the K s filter 

Y 

mago.i 

J 

K s 

50% Completeness level 

Y J K s 

SMC 3_5 

21.9762 

-74.0018 

100% 

21.73 

21.70 

20.62 

22.24 

22.09 

21.02 

SMC 4_5 

21.2959 

-72.9339 

100% 

22.09 

21.72 

20.65 

22.43 

22.15 

21.08 

SMC 5_6 

25.4401 

-71.5879 

92% 

22.02 

21.70 

20.61 

22.52 

22.12 

20.98 

SMC 6_5 

20.4138 

-70.7601 

100% 

22.07 

21.58 

20.59 

22.48 

22.20 

20.93 

BRI 2_3 

33.6941 

-74.0132 

92% 

21.79 

21.62 

20.68 

21.93 

21.71 

20.96 


Table 3. Fundamental parameters of additional SMC clusters. 


Name 

log(t yr 1 ) 

[Fe/H] 

Reference 

BS196 

9.70±0.05 

—1.7±0.1 

i 

HW66 

9.60±0.10 

— 1.3=t0.2 

2 

HW79 

9.70±0.10 

—1.3±0.2 

2 

HW86 

9.15±0.10 

—0.65±0.20 

2, 6 

L100 

9.30±0.10 

—0.7±0.2 

2 

L106 

9.20±0.15 

—0.7±0.2 

2 

L108 

9.30±0.10 

—0.9±0.2 

2, 6 

Llll 

9.15±0.10 

—0.75±0.20 

2, 6 

L112 

9.75±0.10 

—1.1±0.2 

2, 6 

L113 

9.70±0.10 

—1.17±0.12 

3, 4, 6 

L114 

8.15±0.15 

—0.7±0.2 

5 

L115 

8.05±0.10 

—0.7±0.2 

5 


References: (1) |Bica, Santos Sz Schmidt] (|2008|> ; (2) |Piatti et al.| 
(2011); (3) jPiatti et al.j (|2007cj ) ; (4) |Da Costa &: Hatzidimitriou 
(1998); (5) |Piatti et al.| ( |2007b| ; (6) |Parisi et al.| ( |2014[ ). 
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Table 2. Fundamental parameters of the SMC clusters studied. 


Name 

RA 

Dec 

Diameter 

E(B - V) 

log(t yr !) 

Z 

Comments 


(°) 

(°) 

(') 

(mag) 




SMC 3_5 

GKH51* 

022.277 

-73.5635 

0.30 

0.08 



possible non-cluster (1) 

B165 

022.711 

-73.434 

0.45 

0.08 



possible non-cluster 

B166 

022.975 

-73.910 

0.60 

0.08 

8.3 

0.003 


L109 

023.306 

-74.167 

0.90 

0.07 

9.6 

0.003 

9.55±0.10 (2), 9.50d=0.05 (7) 

SGDH-cluster A 

022.308 

-73.532 

0.40 

0.06 



possible non-cluster (1) 

B164 

022.369 

-73.533 

0.60 

0.08 

8.5 

0.003 


GKH54/57* 

022.400 

-73.551 

0.10 

0.08 



possible non-cluster (1) 

GKH29* 

022.400 

-73.558 

0.10 

0.08 



possible non-cluster (1) 

GKH24* 

022.412 

-73.555 

0.05 

0.08 



not resolved (1) 

GKH22 

022.421 

-73.563 

0.40 

0.08 



not resolved (1) 

HW75 

019.372 

-73.570 

0.90 

0.06 

8.3 

0.003 

8.20d=0.40 (6) 

B155 

020.099 

-74.004 

0.50 

0.04 

8.7 

0.003 


SMC 4_5 

HW74 

019.194 

-73.1605 

0.55 

>0.06 

8.8 

0.003 

7.50d=0.20 (6) 

B156 

019.891 

-73.0967 

0.65 

0.06 



possible non-cluster; 7.50d=0.20 (6) 

HW72 

018.923 

-73.167 

0.50 

0.06 

8.7 

0.003 

8.10d=0.20 (6) 

K68,L98 

018.890 

-72.624 

0.70 

0.05 

8.7 

0.003 

8.20d=0.20 (6) 

B162 

020.963 

-73.440 

0.55 

0.06 

7.3 

0.003 


H86-211 

021.191 

-73.425 

0.50 

0.06 



possible non-cluster 

H86-212 

021.313 

-73.501 

0.60 

0.06 



possible non-cluster 

BS282 

021.375 

-73.3895 

0.40 

0.06 



possible non-cluster 

HW77 

020.044 

-72.622 

1.40 

0.06 

9.15 

0.003 


HW78 

020.335 

-73.094 

0.45 

0.06 



possible non-cluster 

HW81 

021.042 

-73.154 

1.00 

>0.3 

7.0 

0.003 


BS176 

021.063 

-73.1685 

0.35 

0.06 

8.3 

0.003 


L110 

023.608 

-72.874 

2.40 

0.06 

9.8 

0.003 

9.80d=0.10 (3), 9.90d=0.10 (7) 

H86-213 

023.673 

-73.2755 

0.40 

0.07 

8.8 

0.003 


HW82 

021.114 

-73.172 

0.70 

0.06 

7.8 

0.003 


HW80 

020.858 

-73.224 

0.70 

0.06 



possible non-cluster 

BS187 

022.754 

-72.851 

0.60 

0.04 

9.3 

0.003 


Skl58 

018.971 

-73.319 

0.70 

0.06 

7.3 

0.003 


Skl57 

018.965 

-73.347 

0.60 

0.06 

7.2 

0.003 


SMC 5_6 

BS188 

023.796 

-71.737 

0.80 

0.07 

9.2 

0.003 


HW84 

025.431 

-71.162 

1.20 

0.07 

9.25 

0.003 

9.20±0.10 (2) 

HW85 

025.615 

-71.2795 

0.70 

0.07 

9.4 

0.003 

9.30±0.10 (2) 

BS189 

025.777 

-71.452 

1.20 

2.80 

7.3 

0.003 


BS190 

025.963 

-71.747 

1.20 

0.50 

8.5 

0.003 


SMC 6_5 

HW73 

019.105 

-71.327 

0.80 

0.05 

8.3 

0.003 

8.15±0.20 (6) 

L95 

018.690 

-71.348 

0.50 

0.05 

8.4 

0.003 

8.30±0.20 (6) 

H86-197 

018.882 

-71.176 

0.50 

0.06 

9.3 

0.006 

9.10±0.10 (5) 

HW67 

018.259 

-70.964 

0.95 

0.06 

9.2 

0.003 

9.20±0.10 (4) 

BS173 

021.017 

-70.327 

1.00 

0.06 



possible non-cluster 

B168 

021.678 

-70.785 

0.50 

0.06 

9.6 

0.003 


IC1708 

021.233 

-71.185 

1.00 

0.07 

9.1 

0.003 

9.10±0.10 (2) 

BRI 2_3 

BS226 

031.425 

-74.381 

0.40 




not recognized 

BS229 

031.952 

-74.4425 

0.50 

0.08 

7.3 

0.003 


BS232 

032.340 

-74.026 

0.50 

0.08 



possible non-cluster 

BS233 

032.658 

-74.156 

1.20 

0.15 

7.3 

0.003 


BS235 

032.962 

-74.119 

0.60 

0.09 



possible non-cluster 

BS239 

033.651 

-73.985 

1.20 

0.11 

8.5 

0.003 


BSBD4 

033.667 

-74.358 

0.60 

0.08 

9.1 

0.003 


BS240 

033.717 

-73.953 

0.35 

0.09 



possible non-cluster 
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Table 2. continued. 


* In Table 3 of|Bica et a l. (2008) the acronym GHK appears instead of GKH (see their table 1). 

References : (1) |Gouliermis, Quanz fc Henning ( 2007} ; (2) |Piatti et al. (|2011| CTj data); (3) Piatti et al.|(|2007c[ CT\ data); (4) |Piatti| 
( |2011b| CTi data); (5) |Piatti Sz Bica ' p012| CT\ data); (6) Glatt, Grebel & Koch ~ p010[ ); (7) |Parisi et al.| ( |2014[ ). 
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Y-Ks 



Y-Ks 



Figure 2. CMDs for stars in the field of HW 77 of tile SMC 4_5: the observed CMD composed of the stars distributed within the cluster 
radius (top left-hand panel); a field CMD for an annulus -outer and inner radii equal to 3.163 and 3.0 times the cluster radius- centred 
on the cluster (top right-hand panel); the cleaned cluster CMD (bottom left). We overplotted boxes for each star in the field CMD to be 
used in the cluster CMD field decontamination (see Section 4 for details). Colour-scaled symbols represent stars that statistically belong 
to the field (P ^ 25%, pink), stars that might belong to either the field or t he cluster (P = 50%, l ight blue), and stars that predominantly 
populate the cluster region (P ^ 75%, dark blue). Three isochrones from pBressan et al.||2012| for log(i yr —1 ), log(£ yr —x ) zb 0.1, and 
the metallicity values listed in Table 2 are also superimposed. The schematic diagram centred on the cluster for a circle of radius three 
times the cluster radius is shown in the bottom right-hand panel. The black circle represents the adopted cluster radius. Symbols are as 
in the bottom left-hand panel, with sizes proportional to the stellar brightnesses. North is up; east is to the left. 
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Figure 3. Spatial distribution of the enlarged SMC cluster sam¬ 
ple (solid circles), where darker solid circles correspond to older 
cluster ages. Catalogued clusters (dots) and ellipses with semima¬ 
jor axes of 2.4° and 6° are overplotted. 











